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Cosmic parallax is the change of angular separation between pair of sources at cosmological 
distances induced by an anisotropic expansion. An accurate astrometric experiment like Gaia could 
observe or put constraints on cosmic parallax. Examples of anisotropic cosmological models are 
Lemaitre-Tolman-Bondi void models for off-center observers (introduced to explain the observed 
ON 1 acceleration without the need for dark energy) and Bianchi metrics. If dark energy has an anisotropic 

equation of state, as suggested recently, then a substantial anisotropy could arise at z < 1 and escape 
the stringent constraints from the cosmic microwave background. In this paper we show that such 
models could be constrained by the Gaia satellite or by an upgraded future mission. 
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I. INTRODUCTION 

In recent times, there has been a resurgent interest towards anisotropic cosmologies, classified in terms of Bianchi 
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solutions to general relativity. This has been mainly motivated by hints of anomalies in the cosmic microwave 
background (CMB) distribution observed on the full sky by the WMAP satellite [H, 0, H 3 • While the CMB is very 
well described as a highly isotropic (in a statistical sense) Gaussian random field, recent analyses have shown that 
local deviations from Gaussianity in some directions (the so called cold spots, see Q) cannot be excluded at high 
confidence levels. Furthermore, the CMB angular power spectrum extracted from the WMAP maps has a quadrupole 
power which appears significantly lower than expected from the best-fit cosmological model Q. Several explanations 
for this anomaly have been proposed (see e.g. @, 0, B B including the fact that the universe is expanding with 
different velocities along different directions. While deviations from homogeneity and isotropy are constrained to be 
very small from cosmological observations, these usually assume the non-existence of anisotropic sources in the late 
universe. Conversely, as suggested in [ll|, [lj, [13, [lj, [l5| , dark energy with anisotropic pressure acts as a late-time 
source of anisotropy. Even if one considers no anisotropic pressure fields, small departures from isotropy cannot be 
excluded, and it is interesting to devise possible strategies to detect them. 

The effect of assuming an anisotropic cosmological model on the CMB pattern has been studied by [3, E3, 0, 
IHi I20L [2l| . The Bianchi solutions describing the anisotropic line element were treated as small perturbations to 
a Friedmann-Robertson- Walker (FRW) background. Such early studies did not consider the possible presence of a 
non-null cosmological constant or dark energy and were upgraded recently by [13, [23| . 

One difficulty of the anisotropic models that have been shown to fit the large-scale CMB pattern is that they have 
to be produced according to very unrealistic choices of the cosmological parameters. For example, the Bianchi Vllh 
template used in [23] requires an open universe, an hypothesis which is excluded by most cosmological observations. 
An additional problem is that an inflationary phase - required to explain a number of feature of the cosmological model 
- isotropizes the universe very efficiently, leaving a residual anisotropy that is negligible for any practical application. 
These difficulties vanish if an anisotropic expansion takes place only well after the decoupling between matter and 
radiation, for example at the time of dark energy domination [ill, El, E, fl^-flH|. 

The effect of cosmic parallax [131 has been recently proposed as a tool to assess the presence of an anisotropic 
expansion of the universe. It is essentially the change in angular separation in the sky between far-off sources, due 
to an anisotropic expansion. This all-sky change in separations can be used as a tracer of anisotropic behaviour of 
the spacetime metric. This effect has been investigated in the context of Lemaitre-Tolman-Bondi (LTB) models with 
off-center observers In this paper we study the cosmic parallax in Bianchi I metrics (see also [25] for ellipsoidal 
universes). We will show that, since the cosmic parallax traces the geodesic of the metric to the present time, it can 
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be used to constrain the late anisotropic behaviour induced, for example, by the above mentioned anisotropic dark 
energy models. This makes it a valuable tool with respect to primary CMB anisotropies which are frozen at z ~ 1000. 

While finalizing this paper another work analysing the cosmic parallax in Bianchi I models appeared [25]. We will 
discuss the main differences with this work later on. 



II. COSMIC PARALLAX IN BIANCHI I 



We consider a class of homogeneous and anisotropic models where the line element is of the Bianchi I type, 

ds 2 = -dt 2 + a 2 {t)dx 2 + b 2 {t)dy 2 + c 2 {t)dz 2 . (1) 

The expansion rates in the three Cartesian directions x, y and z are defined as Hx = a/a, Hy = b/b and Hz = c/c, 
where the dot denotes the derivative with respect to coordinate time. In these models they differ from each other, 
but in the limit of Hx = Hy = Hz the flat FRW isotropic expansion is recovered. Among the Bianchi classification 
models the type I exhibits flat geometry and no overall vorticity; conversely, shear components Tix,y,z = Hx,y,z / H — 1 
are naturally generated, where H is the expansion rate of the average scale factor, related to the volume expansion 
as H = A/A with A = (abc) 1 '*. 

Cosmic parallax is the temporal change of angular separations between distant sources in the sky caused by large 
scale anisotropic expansion |24l |. The sources are assumed to trace the evolution of the cosmic expansion (see for 
example also [26j]); since the parallax induced by peculiar velocity is randomly distributed, it can be averaged out of 
a large sample and, in addition, decreases with distance from the observer [24 |. 

For an off-centre observer in a LTB model the cosmic parallax is a pure dipole signal in the sky that might be affected 
by systematic noises like the observer's peculiar velocity and acceleration (the latter induces aberration changes), even 
though an observational strategy using different source samples at different redshifts (say, within and outside the void) 
would help to discriminate between them. In homogeneous and anisotropic models like Bianchi I we expect the signal 
to have a different angular distribution, hence being even more predictive. 

Let us consider two sources A and B in the sky located at physical distance from us observers 0[a,b] = 
(X, Y, Z)yA.B\ = (R sin 8 cos (j), Rsinds'mcj), R cos 6)[a,b] , where R = V 1 X 2 + Y 2 + Z 2 and {8,4>) are spherical angular 
coordinates. Their angular separation on the celestial sphere reads 

OA - OB = cos 7 = cos 8 a cos 6b + sin 6a sin 6 b cos A<fi, (2) 

with A(j) = {4>a — 0s). From now on we will mark spatial separations and temporal variations with A and A t symbols, 
respectively. If A t j ^ then a cosmic parallax arises. In a homogeneous and isotropic model (like FRW) the geodesic 
are radial and sources are subjected to the same radial expansion rate, keeping their angular separation constant. 
On the other hand, if the expansion is anisotropic their spherical coordinates change dissimilarly in time leading to a 
modification of their angular separation: 

— sin7A t 7 = sin 6 a cos 6B(A t 6B cos Acj) — A t 6 a) + cos 8 a sin 6b {A t 6 a cos A(f) — A t &B) (3) 

+ sin 8 a sin 8b sin Aifi(At4>B — A t 4> a)- 

In the limit A t (j)A = A t (j)B = 4>A = <i>B = the relative motion is constrained on the (X,Z) plane and the cosmic 
parallax reduces to (At 0,4 — A t 8s) (see Fig. [I]). Similarly, on the (X,Y) plane the signal is {A t <j)A — A t 4>B) (see Fig. [2]). 
Both in Fig. [TJ and [2] we allowed the shear parameters at present to appreciably deviate from 0. This explains why the 
cosmic parallax is few orders of magnitude larger than the one in [23]. The main motivation for this will be presented 
in the first paragraph of Section IIVI 

The signal is in general a combination of both the anisotropic expansion of the sources themselves and the change 
in curvature induced by the shear on the photon path from the emission to the observer. In inhomogeneous and/or 
anisotropic models photons follow trajectories that, in general, are not radial. However, while in LTB models this 
effect on the cosmic parallax is enhanced by inhomogeneity (although a FRW description of null geodesic has been 
shown to be fairly good approximation [24j), in Bianchi I models we consider in this paper the geodesic bending 
for a single source amounts at most to about 7% (see Appendix A), which allow us to adopt the straight geodesies 
approximation. 

The spherical angular coordinates are related to the Cartesian coordinates via 4> — arctan(F/A) and 8 — 
arccos (Z/^f X 2 + Y 2 + Z 2 ). Therefore their time evolution can be written as 
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Figure 1: Cosmic parallax in Bianchi I models as a function of 9 for <j> = A(f> = and AO — 90° (this setting corresponds to the 
(X,Z) plane), h x = 0.71, h y = 0.725, h z = 0.72, i.e. E x = -0.012 and E y = 0.009). The time interval is AT = lOyrs. 




Figure 2: Cosmic parallax in Bianchi I models as a function of <f> for = 90°, Ad = and A<p = 90° (this setting corresponds 
to the (X,Y) plane), h x = 0.71, h y = 0.725, h z = 0.72, i.e. E x = -0.012 and E y = 0.009). The time interval is AT = lOyrs. 

where we made use of the Hubble law in the three cartesian directions A t X, Y, Z = X,Y, Z ■ H$x,Y,zAt, valid for 
small time span At (which we assume of the order of decades) relatively to the cosmic time. 
In spherical coordinates Eqs. <j4jt5j) can be expressed as 

A t (f> = ^^(Sor-Sox)-ffoAt (6) 



A t 6 = ^ 



(3(£ x + Z QY )+ cos 2<£(£ x - Soy )J H At, (7) 

where 1Box,Y,z are the shear components at present as defined at the beginning of this section, satisfying the transverse 
condition Sox + Soy + Soz = 0. 

Equations |(6][7j) describe a pure quadrupole signal in the 4> an d 9 coordinate, respectively. This functional form of 
the signal is exactly the same as the one expected for the first non-vanishing multipole expansion of the CMB large 
scale relative temperature anisotropies in Bianchi I model [l8| (remember we are neglecting all peculiar velocities, 
including our own). By combining them into the full spherical distance formula ([3]) the resulting cosmic parallax 
signal obviously exhibits a more complicated shape depending on the location on the sky. 

Considering two sources with an initial angular separation such that (#,</>) b = {9,4>)a + A(0, 0) and substi- 
tuting Eqs. (JGKTJ) in Eq. ((3j) we gain the full expression for the cosmic parallax in spherical coordinates A t 7 = 
Atj(6, <j), A6, A(j>, Y.ox , S y , H , At), where the only further conjecture is that H does not vary in At. At first order, 
this seems reasonable for the time intervals under consideration. Notice that the signal turns out to be independent on 
the redshift: source pairs along the same line of sight undergo the same temporal change in their angular separation. 
This means that aligned quasars would stay aligned. In general of course the number density of quasars will change 
so that the number counts should show some level of anisotropy. This could provide an additional constraint on 
anisotropic expansions: we will discuss briefly this possibility in Sect. V. 

If there were no anisotropies present at last scattering of course a late time anisotropic expansion would marginally 
affect the CMB via a late time direction dependent integrated Sachs- Wolf. Notice that the positional shift of the 
sources themselves is a completely different signal with respect to the bending of light ray during propagation time. 

In Fig. OS, Mollweide projections on the sky of the cosmic parallax signal with respect to a fixed source located on 
the north pole and on the (X,Y) plane are shown. As expected, when the source is at an equatorial position the 
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Experiment N s 



G acc 



At 



Gaia 



500,000 50^as 5yrs 



Gaia+ 



1,000,000 5Aias lOyrs 



Table I: Specifications adopted for Gaia-like and Gaia+ experiments, where N s is the total number of sources, o acc is the 
experimental astrometric accuracy and At is the time interval between two measurements. 

symmetry with respect to the (X,Y) plane is preserved, while when the source is at the north pole a symmetry with 
respect to the (X,Z) plane emerges. In a FRW universe the components of the shear simultaneously vanish and so 
does the cosmic parallax. 



As a next step, we would like to give an insight on whether accurate future satellite astrometry mission will be 
able to put constraints on the anisotropy parameters that are competitive with CMB quadrupole constraints [l§1 |. 
An astrometry mission like Gaia will detect around 500,000 quasars in its 5 years flight time with positional error 
Oacc =10-200 fi&s [13, H3|- Attributing to a Gaia-like experiment the capability of detecting the quasar angular 
positions at two different time separated by At w lOyrs (i.e. conceiving the possibility of two separated missions or 
just a longer one) we can adopt its instrumental characteristics to perform a Fisher matrix analysis. In these Bianchi 
I models the cosmic parallax signal depends on four parameters: the average Hubble function at present, the time 
span and the two Hubble normalized anisotropy parameters at present. However, for the allowed range of values, 
contours in the (Sox, Soy) frame do not depend on the value of Hq, Stretching the time interval between the two 
measurements or improving the instrumental accuracy would instead have an impact on the final constraints. In order 
to analyse these dependencies we make use of the Fisher formalism, namely the Fisher Matrix defined as 



where all separations are taken with respect to a reference source and index I runs from 2 up to the number of quasars 
N s to take into account the spherical distances to all other sources. In fact, one should notice that the Gaia accuracy 
positional errors are obtained having already averaged over 2N S coordinates. 

We simulated a catalogue of up to 1,000,000 quasars with angular positions (6, <fi) randomly generated from a 
uniform distribution on the celestial sphere. We then used the covariance matrix Cy = (-Fy) -1 to construct the error 
ellipses with 1 and 2a contours in the (Sox, Soy) plane. The positional accuracies should have a mild dependence on 
the magnitude. The quasars are expected to have magnitudes ranging from 12 to 20 and, correspondingly, accuracies 
from 10 down to 200 ^tas, as pointed out in [28]. We could have weighted our non-redshift dependent signal with 
accuracies that are function of magnitude. However, for simplicity we adopt a single representative average accuracy 
of a a cc =50 /ias; it is immediate to rescale the final errors to a different accuracy. We also perform the calculation for 
an enhanced Gaia-like mission dubbed as Gaia+ (see specifications in Table ITTT|) . 

The Fisher error ellipse are shown in Fig |4j the constraints turn out to be of the same order of magnitude of the 
CMB limits on the shear at decoupling. The la errors on S x and S y are 8.3 • 10 -4 and 6 • 10~ 5 for Gaia and Gaia+, 
respectively. Although our null hypothesis was chosen to be the friedmannian isotropic expansion (S x = S y = 
and Hq — 72 km/s/Mpc), due to the linear dependence of the signal on the shear parameters, a change of the fiducial 
model corresponds to a simple translation of the same ellipse in the frame. Gaia data processing is incredibly complex 
and the experimental covariance matrix will probably be at the end non-diagonal. However, the Gaia collaboration 
have not provided yet a quantification of these correlations and giving a formal status to it is beyond the scope of 
this paper. 



III. COSMIC PARALLAX FORECASTINGS 
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Figure 3: Mollweide contour plot for cosmic parallax in Bianchi I models for one source fixed at two different location in the sky. 
Upper panels show the signal for ellipsoidal models (hoz = 0.72 and hox = hoy = 0.71), while in lower panels hoy = 0.725. 
Lighter colours correspond to higher signal and on the horizontal and vertical axes angular coordinates vary in the range 
4> : [0, 2ir] and 6 : [0, it], respectively. The time interval is AT = lOyrs. 




Figure 4: Fisher contours of Cosmic parallax for Gaia and Gaia+ specifications (dashed and solid lines, respectively). The 
double contours identify la and 2a regions for AT = lOyrs. 



IV. COSMIC PARALLAX INDUCED BY DARK ENERGY: AN EXAMPLE 

The CMB provides very tight constraints on Bianchi models at the time of recombination H3, HH of order of 
the quadrupole value, i.e. ~ 10~ 5 . Usually, in standard cosmologies with a cosmological constant the anisotropy 
parameters scale as the inverse of the comoving volume. This implies an isotropization of the expansion from the 
recombination up to the present, leading to the typically derived constraints on the shear today, namely ~ 10~ 9 -^10~ 10 
(resulting in a cosmic parallax signal of order 10 _4 ^as). However, this is only true if the aforementioned parameters 
are monotonically decreasing functions of time, that is if the anisotropic expansion is not generated by any anisotropic 
source arising after decoupling, e.g. vector fields representing anisotropic dark energy [la ]. 

Motivated by this, we apply the cosmic parallax to a specific anisotropic phenomenological dark energy model 
in the framework of Bianchi I models [TJ E2] (we refer to these papers for details) . This description allows semi- 
analytical calculations and represents in a fairly conservative approach more complicated anisotropic models. Here 
the anisotropic expansion is caused by the anisotropically stressed dark energy fluid whenever its energy density 
contributes to the global energy budget. If the major contributions to the overall budget come from matter and dark 
energy, as after recombination, their energy-momentum tensor can be parametrized as: 

T ( m ) v = Aiag(-1, w m ,w m ,w m )p m (9) 
r ( D E)i/ = diag(-l, w + 3(5, w + 3-f)p DE , (10) 
respectively, where w m and w are the equation of state parameters of matter and dark energy and the skewness 
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parameters 6 and 7 can be interpreted as the difference of pressure along the x and y and z axis. Note that the 
energy- momentum tensor (fTUj) is the most general one compatible with the metric HJ [12]. Two quantities are 
introduced to define the degree of anisotropic expansion: 

R = (a/a - b/b)/H = S x - E„ , 

S = (a/a-c/c)/H = 2E X + £ a . [ ' 

The reason why the cosmic parallax is allowed to be few orders of magnitude larger than the one in [25| is based on 
the presence of this anisotropic source arising after decoupling. In particular, the value 6 = —0.1 is not completely 
excluded by supernovae data, since it lies on the 2a contours of the gamma-delta plane, if a prior on w and f2 m is 
assumed [12 [. More phantom equation of state parameters and/or larger matter densities allow for larger value of 
delta. In addition, and more in general, time dependent delta and gamma functions, mimicking for example specific 
minimally coupled vector field with double power law potential, can escape these constraints. Our purpose is to use 
this parameterization to model a very late-time evolution of the shear, which is the reason why we linearised the 
dynamical solutions around the critical points as denoted in the following paragraphs. 

Considering the generalized Friedmann equation, the continuity equations for matter and dark energy and no 
coupling between the two fluids, the derived autonomous system reads 1 (ill. fl2|: 

U' =U(U - 1)[ 7 (3 + R-2S) + 6(3 -2R + S) + 3(w - w m )} 

S' =i(9 - R 2 + RS - S 2 ){S[U(6 + -/ + W- w m ) + w m - 1] - 67E/} (12 ) 

R' =i(9 -R 2 + RS- S 2 ){R[U(6 + 7 + w - w m ) + w m - 1] - 6 SU}, 

where U = pde/(pde + Pm) and the derivatives are taken with respect to log(A)/3. In what follows we will consider 
for simplicity that w m = 0, i.e. pressureless matter . System lfl2|) exhibits many different fixed points, defined as the 
solutions of the system S' = R' = U' = 0. Beside the Einstein-de Sitter case (R* — S* = U* = 0), the most physically 
interesting for our purposes are the dark energy dominated solution 

R. = - x r, S.= - s 67 - , U* = l, (13) 

(5 + 7 + w-l J + 7 + W-I v ' 

and the scaling solution 

= 3(5((5 + 7 + w) = 37((5 + 7 + w) = ^ + 7 + ^ 

2(J 2 -5 7 + 7 2 )' * 2(6 2 - 6 1 + 7 2 ) ' * w 2 - 3(7 - 5) 2 + 2w(j + 6) ' 1 ' 

in which pw/p m — const., i.e., the fractional dark energy contribution to the total energy density is constant. The 
latter is positive if the numerator and the denominator in the expression for U* are either both positive or both 
negative; moreover we should ensure the condition U* < 1, If the numerator is positive then w > —(j+5) is required; 
the denominator is then positive only if w > Wi = — (7 + 6) + ^/( 7 + 5) 2 + 3(7 — 6) 2 , with wi being positive. Hence 
the most interesting case is when they both are negative, which translates into these conditions: 1) w < — (7 + 6); 
2) w > W2 = — (7 + 5) — v(7 + 5) 2 + 3(7 — S) 2 , where now W2 is negative. There are other conditions one should 
impose to such solutions, in particular stability: we refer the reader to Ref. [ill, EH for the details. All the specific 
cases discussed below satisfy these conditions. For simplicity we also assume w = — 1 and 7 = 0. 

The cosmic parallax constrains the anisotropy at present, when the dark energy density is of order 74%, hence not 
yet in the final dark energy dominant attractor phase 1(13)1 . Therefore it must be either on its way to such a stage 
or, alternatively, on the scaling critical solution lfT4)l . We discuss separately these two alternative scenarios. In the 
scaling case, in order not to produce a too long accelerated epoch in the past, we ensure that we just entering the 
accelerated regime. 

We map our Fisher matrix $8$ with the same experimental specifications as in Table ITTT1 into the new parameter 
space p = (6, 7). Our intention is to infer the order of magnitude of the constraints this new cosmological tool 
would be able to put on the dark energy skewness parameters. We apply the parameter transformation F' — A T FA, 
where Ay = dT, i/dpj. If we assume that the system has just entered the scaling solution, the critical point lfl4)l 
approximately describes the dependence of the dark energy anisotropy on the skewness parameters at present. In this 



1 Notice that in [IJ there is a spurious factor x in the phase-space equations (8). 
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Figure 5: Projected Fisher contours for the skewness dark energy parameters for Gaia (upper panels) and Gaia+ specifications 
(lower panels). The double contours identify la and 2a regions for AT = lOyrs. The dashed lines represent the case of an 
ellipsoidal universe with w = — 1, Uo = 0.74 and S = —0.1 (Ro =■ 0.2) approaching the dark energy dominated critical point 
(where (7 = 1 and R* ~ 0.3), while the solid lines represent an ellipsoidal universe that has just entered the scaling regime, 
with w = -1, Uo =0.74 and 5 = 0.5 (R* ~ -0.5). 




-3 -2-10 1 

Log(abc) 

Figure 6: Cosmological evolution of R and U for the two cases considered (both with w = — 1 and Uo = 0.74): {5 = —0.1 , Ro — 
0.2} (solid lines) and {8 = 0.5 , Ro ~ -0.5} (dashed lines). 

case we choose for the fiducial model 5 = 0.5 and 7 = 0, namely an ellipsoidal Universe with Rq = —0.5 (see Fig [6]) 
and we choose initial conditions such that Uo — 0.74. The final attractor value is as expected rather close, Uq — 0.67. 
The error contours are shown in the right panels of Fig. [H for both Gaia and Gaia+ configurations. 

Conversely, if the expansion is driven towards a future dark energy dominated solution, equations (fl~3|) do not 
represent the anisotropy parameters at present (see Fig [6]) . In order to derive a more appropriate functional form for 
them, we solved the linearized system lfl2|) around solution (jT3j) and we fixed \og(A) = to select the present values. 
For this second case, results are shown in the left panels of Fig. [H Notice that Fig [6] depicts a late time expansion 
history and aims just at illustrating the trend towards the critical points from different values of the anisotropy 
parameters (namely R — 3 in an earlier sta ge, but it might set to be vanishing at decoupling by time dependent 
skewness parameters in specific models - see [12l|). 



8 



The constraints on 7 and 5 are in the range 10~ 3 -j- 10 -4 . The current limits from SNIa data are then 2 to 3 
orders of magnitude weaker [12] and, even if the number of supernovae will substantially increase in the near future, 
it might be hard to improve the constraints at such a high level because of the integral dependence of the luminosity 
distance on the skewness parameters. Therefore the cosmic parallax seems to be an ideal candidate for testing the 
anisotropically stressed dark energy. 

The forecastings we presented in this section do not include possible systematic effects. In our Fisher analysis 
we just took into account the statistical errors. Several spurious effect must be considered by the time real data is 
available. For instance, the peculiar velocity of the objects need to be considered, although averaging it over a large 
sample of uncorrelated objects it should be possible to eliminate such a form of bias. Furthermore, this effect decreases 
with increasing angular diameter distance to the object [H[. The main source of noise however could be due to the 
aberration change induced by our own motion. Fortunately, both the aberration change and the observer peculiar 
velocity signal have a dipolar signature whereas cosmic parallax from Bianchi I models results in a superposition of 
quadrupoles. Other minor effects like the temporal changes of local Jensing and microlensing (a parallax disturbances 
of few nanoarcseconds are expected due to the weak microlensing [2s|). were not taken into account, as they require 
a more detail analysis, beyond the scope of this paper. 

V. DISCUSSIONS 

Any anisotropy will leave an imprint on the angular distribution of objects that are able to trace cosmic expansion. 
If such an anisotropy is present before the last scattering surface, the CMB map will also be affected. The temperature 
field will carry extra anisotropics mainly caused by the angular dependence of the redshift at decoupling. By resolving 
geodesic equations and expanding temperature anisotropics in spherical harmonics, it is straightforward to relate the 
low multipole components to the eccentricities of the model |12}. In Bianchi I models the first notable multipoles 
related to the CMB are the monopole and the quadrupole. The observed value of the latter puts constraints on the 
shear at last scatteringoi order 10 -5 , taking into account the cosmic variance. These constraints can be mapped into 
either magnetic field [9( or anisotropic dark energy limits, depending on the source that gives rise to the anisotropy. 
In addition, one expects the eccentricities to be non-vanishing if the expansion has been somewhat anisotropic at 
decoupling. However it is in principle also possible to escape detection from CMB if each scale factor has expanded 
the same amount since last scattering, no matter how anisotropically. Nonetheless, in all these analysis the anisotropy 
pattern is directly added to the intrinsic standard FRW perturbations, a simplistic way of treating the signal at 
large scale. In this still exploratory stage of analysis it is worth stressing that CMB is indeed a very powerful 
constraint on the shear at the time of decoupling, but with almost no direct impact on late time expansion history. 
Complementary to that, cosmic parallax, namely the temporal change of angular separation of distant sources, is a 
direct and potentially powerful test of anisotropy at small redshifts and at present. 

The anisotropic stress of dark energy is expected to have a leading role in the generation of anisotropy at late times. 
It can be parameterized by skewness parameters in the stress-energy tensor formulations, which may be constant 
or time dependent functions. For example a minimally coupled vector field satisfying quadrupole constraints was 
presented in [12] • The only way to test these models is to use either the angular dependence of the magnitude or the 
angular distribution of objects in the sky at recent time, i.e. either distant source angular distribution or the real-time 
cosmic parallax, the two relying on different techniques and having independent systematics that complement each 
other. 

We adopted a simple phenomenological model to describe late-time expansion of the universe, filled by pressureless 
matter and an anisotropically stressed dark energy component with two extra degrees of freedom, namely the skewness 
parameters. No matter what the evolution of both energy densities and shears was at high redshifts, we have shown 
that Gaia will be able to constrain the skewness parameters up to 10~ 3 -j- 10~ 4 at 2cr, comparable to CMB tests at 
decoupling time, and 2-^3 orders of magnitude better than current supernovae la limits [12| (a Gaia+ experiment 
would improve them by one order of magnitude) . 

In this paper we discussed in detail the real-time technique; before concluding we comment briefly on the possibility 
of testing anisotropy through the accumulated effect on distant source (galaxies, quasars, supernovae) distribution. 

If sources shifts by as much as 0.1/ias/year during the dark energy dominated regime, then the accumulated shift 
will be of the order of 1 arcmin in 10 9 years and up to fraction of a degree in the time from the beginning of acceleration 
to now. If the initial distribution is isotropic, this implies that sources in one direction will be denser than in an 
orthogonal direction by roughly 1/90 » 10 -2 . This anisotropy could be seen as a large-scale feature on the angular 
correlation function of distant sources, where we expect any intrinsic correlation to be negligible. The Poisson noise 
become negligible for N ^> 10 4 : for instance, a million quasars could be sufficient to detect the signal. Although 
the impact of the selection procedure and galactic extinction is uncertain, this back-of-the-envelope calculation shows 



9 



that the real-time effect could be complemented by standard large-scale angular correlation methods 2 . 

While finalizing this manuscript another work analysing the cosmic parallax in Bianchi I models came out [HI- 
Our work differs in many aspects. In [HI the authors focused on the shear in models that isotropize, that is on 
solutions of the shear dynamical equations that are decreasing function of time. Hence the fact that they find a 
signal substantially lower than ours is not surprising. In fact, in [25| the background is described by a ACDM model 
where the non-FRW quantities are driven by a constant equation of state. The analysis is restricted to ellipsoidal 
universes, where the dependence on the azimuthal angle is dropped and the signal is a pure quadrupole. Furthermore 
in order to accomplish forecasting we have performed a Fisher Matrix analysis of the signal contemplating two different 
experimental sets, Gaia and Gaia+. 

More in general, we have shown that, differently from LTB models with off-centre observers [Hj], the cosmic parallax 
signal in Bianchi I models is a combination of two quadrupole functions of the two angular coordinates. Since the 
most important systematic noises, caused by peculiar velocities and aberration changes, have a dipolar functional 
form, Bianchi I models seem to be ideally testable, though even in LTB models specific observational strategy aiming 
at distinguishing the signal from the noise are possible [241 ]. 

Assuming that null geodesies are radial, we have provided an analytical expression for the cosmic parallax in 
general Bianchi I models. This assumption is motivated by a direct numerical calculation of the geodesic for a source 
at redshift z — 1 (see Appendix A). 

CMB and cosmic parallax detect anisotropy at two different times and, from an observational point of view, are 
completely independent on each other: combining them together one will have the opportunity to reconstruct the 
evolution of the anisotropy and test with high accuracy the Copernican Principle. 
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Appendix A: GEODESIC EQUATIONS 



Restricting for simplicity to two dimensions, particularly to the (X,Y) plane (where — 7r/2), we now want to 
check whether neglecting the curvature of null geodesic equations considerably affects our results. Photons follow 
trajectories that are described by the ensuing equations: 

X" = -2H x t'p (Al) 
Y" = -2H Y t'q (A2) 

z' = - { -^^(a 2 H xP + b 2 H Yq ), (A3) 

with the additional constraint t' 2 = a 2 p + b 2 q (note that here z is the redshift, not to be confused with the coordinate 
Z). Here primes denote derivative with respect to the affine parameter A. 

In order to solve system (|AlHA3jl , the scale factors as functions of time are required and hence one has to couple 
to it the dynamical equations. Since we are integrating backward from the observer position to the source location, 
typically at redshift of order 1, we need to evaluate these functions in these redshift range. We adopt the linearized 
solution of the dynamical system (fl2|) in the vicinity of the critical point (fl3|l : in this way we take into account the 
effect of the shear arising from an anisotropically stressed dark energy. The linearized equation for the anisotropy 
parameters with (5 = 7= —0.1 and w = — 1 are: 

S(A) = 0.27-0.4(L4- 3 - 50 + 0.38/r 3 - 27 (A4) 
R[A) = 0.27- 0.07A" 3 ' 27 (A5) 
U(A) = 1 - 0.26A" 3 - 27 . 

The full set of equations with S(A), R(A) and U(A) is then : 



2 We are indebted to an anonymous referee for this suggestion. 
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Figure 7: Geodesic bending for a source at z = 1 in the "X-Y" plane for the the dynamical system lfl2|) in the vicinity of the 
critical point l113[l in four different cases: photons arriving at (pa = 90° (along the F-axis), <J>o = 45°, (f>o = 22.5° and <f>o = 0° 
(along the X-axis). The total deviation is always less than 7%, which validates the "straight geodesies" approximation. 



(j) 2 [l-l(S 2 + R 2 -RS)] = H$(Q m0 A- 3 + Q deA-^ 1+w+s ^). 
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Results for a source located at z=l are shown in Fig. [7l 
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